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Abstract 

The Pierre Auger Observatory for cosmic rays provides a laboratory for studying fundamental 
interactions at energies well beyond those available at colliders. In addition to hadrons or 
photons, Auger is sensitive to ultra-high energy neutrinos in the cosmic radiation and models 
for new physics can be explored by observing neutrino interactions at center-of-mass energies 
beyond the TeV scale. By comparing the rate for quasi-horizontal, deeply penetrating air 
showers triggered by all types of neutrinos with the rate for slightly upgoing showers generated 
by Earth-skimming tau neutrinos, any deviation of the neutrino-nucleon cross-section from the 
Standard Model expectation can be constrained. We show that this can test models of low- 
scale quantum gravity (including processes such as Kaluza-Klein graviton exchange, microscopic 
black hole production and string resonances), as well as non-perturbative electroweak instanton 
mediated processes. Moreover, the observed ratios of neutrino flavors would severely constrain 
the possibility of neutrino decay. 



I. INTRODUCTION 



The Pierre Auger Observatory is the largest cosmic ray detector in the world p|. 
Currently under construction at Malargiie (Argentina), it has begun taking data and 
already accumulated an exposure comparable to previous experiments such as AGASA 
and HiRes 0. 

In addition to studying the highest energy cosmic rays, Auger is also capable of ob- 
serving ultra-high energy cosmic neutrinos At present, the AMANDA telescope at 
the South Pole holds the record for the most energetic neutrino interactions observed 
these events have energies up to ~ 10^ GeV and are consistent with the predicted spec- 
trum of atmospheric neutrinos. Auger, by contrast, is expected to detect neutrinos with 
energies above ~ 10^ GeV. The ability to study neutrino interactions at such high ener- 
gies will open a unique window on possible physics beyond the Standard Model (SM) of 
strong and electroweak interactions. 

A variety of models have been proposed in which neutrino interactions become sub- 
stantially modified at very high energies, the most interesting being models of low scale 
quantum gravity (involving the exchange of Kaluza-Klein (KK) gravitons, production 
of microscopic black holes, and excitation of TeV-scale string resonances), and models 
featuring non-perturbative electro-weak instanton induced interactions. 

Moreover the neutrino flavor ratios predicted by the standard oscillation phenomenol- 
ogy can be modified in propagation over cosmological distances if processes such as neu- 
trino decay occur. Auger is expected to detect the 'cosmogenic' neutrino flux from 
interactions of extragalactic ultra-high energy cosmic rays with the cosmic microwave 
background. Thus, it will be sensitive to such effects, being capable of measuring the 
flux of ultra-high energy tau neutrinos in addition to the overall neutrino flux. 

In this article, we explore the relevant phenomenology and quantify the sensitivity of 
Auger to such new physics.^ In Secim we describe the Auger experiment and its ability to 
detect quasi-horizontal neutrino-induced air showers, as well as up-going showers induced 
by Earth-skimming tau neutrinos. In Sec. Illll we discuss possible sources of cosmic ultra- 
high energy neutrinos. In section IIV| we infer the sensitivity of Auger to the neutrino- 
nucleon interaction cross-section and to the flavour content of the ultra-high energy cosmic 
neutrino flux. In Sec. |3we consider speciflc models of physics beyond the SM and their 
signatures in Auger. We present our conclusions in Sec. IVII 



II. ULTRA-HIGH ENERGY NEUTRINOS AT AUGER 
A. The Pierre Auger Observatory 

Auger is a hybrid ultra-high energy cosmic ray detector, with a ground array of water 
Cerenkov detectors sampling air shower particles, overlooked by air fluorescence detector 
telescopes which observe the longitudinal development of the showers When com- 
pleted in 2005-06, the Southern hemisphere site will have 1600 detectors on the ground 
covering 3000 km^, and 4 fluorescence telescopes having 6 detectors each. A similar 
facility has been proposed for a Northern hemisphere site in Colorado (USA). In our 
calculations we consider only a single site and focus on the ground array. 



For a review of exotic neutrino interactions and their signatures in high-energy cosmic neutrino tele- 
scopes such as IceCube, see Ref. H. 



B. Quasi-Horizontal, Deeply Penetrating Showers 



At sufficiently high energies cosmic neutrinos can trigger atmospheric air showers 
similar to those due to high energy cosmic rays (hadrons or photons). However, unlike 
ordinary cosmic ray showers which are initiated near the top of the atmosphere, those 
generated by neutrinos can be initiated at any depth since the interaction cross-section 
is much smaller, hence the probability of interaction per unit length is approximately 
constant. Neutrino induced showers can thus be distinguished from cosmic ray showers 
by requiring that they be deeply penetrating. This is most useful for distinguishing the 
two kinds of showers, because within ~ 20° of the horizon the electromagnetic component 
of hadron- induced showers is completely absorbed before reaching the detector. 
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FIG. 1: Angular distribution of neutrino- induced showers expected to be observed by the Auger 
ground array with different selection criteria. The upper solid line indicates all showers gen- 
erated by cosmic neutrinos over the range of zenith angles shown, while the upper and lower 
dashed lines correspond to the cases when the shower is initiated at a depth exceeding 1000 
and 2000 g/cm^, respectively. The lower solid line corresponds to the case when the shower is 
initiated at a depth exceeding 2000 g/cm^ and within 2000 g/cm^ from the detector — this last 
set of cuts is adopted in our calculations. For the four cases shown, the fraction of events which 
survive the various cuts are 1.0, 0.80, 0.60 and 0.33, respectively. We have adopted a neutrino 
spectrum oc E^"^ , saturating the Waxman-Bahcall flux bound, see Eq. 

The rate of neutrino-induced showers expected to be observed in an experiment such 
as Auger can be written as 
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where ATobs is the observation time, Aa is Avogadro's number, dauN/dy is the differ- 
ential neutrino-nucleon cross-section, y is the inelasticity, 6^ is the zenith angle, A± is 



the cross sectional area of the experiment as seen from a given zenith angle, X is the 
atmospheric depth of the interaction (the atmospheric mass per unit area), Egh is the 
total energy dissipated in the shower, and dNi,/dE,^ is the incoming neutrino flux. The 
function P[Esh,Oz,X] is the probability of the experiment detecting a shower created at 
an atmospheric depth X of energy E'sh, at a zenith angle O^- In order to ensure that 
such a shower can be distinguished from one initiated by a hadron or photon primary, 
we require that for the Auger ground array: 

• The zenith angle, 6^ > 70°, 

• The neutrino penetrates at least 2000 g/cm^ into the atmosphere before interacting, 

• The interaction takes place within 2000 g/cm^ of the detector. 

This third requirement is included due to the difficulty in reconstructing events beyond 
this range at Auger. At 6z = 70°, the total path length traversed before the shower hits 
the Earth's surface is Xground ~ 3000 g/cm^, hence relatively little atmosphere is present 
in which a neutrino primary can interact and be distinguished from an ordinary cosmic 
ray shower. At zenith angles larger than 85° however, the slant depth exceeds 10000 
g/cm^. In Fig. Q we show how these selection criteria will affect the observed angular 
distribution of quasi-horizontal showers in Auger. 
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FIG. 2: Spectra of quasi-horizontal, deeply penetrating, neutrino induced showers as would 
be seen by Auger with the cuts shown in Fig^ The solid and dashed lines correspond to 
the cosmogenic flux and the Waxman-Bahcall flux (see section IIIII for details) which yield, 
respectively, 0.07 and and 0.22 events per year. For the latter case, the 'bumps' at ~ 10^'^ GeV 
and ~ lO*^'^ GeV correspond to NC and CC interactions respectively, while the W~ resonance 
can be seen at 6.3 x 10^ GeV. 

If the neutrino-nucleon cross-section is enhanced well above the SM prediction at ultra- 
high energies, it is possible that this depth of atmosphere will significantly attenuate the 
cosmic neutrino flux; to account for this, an additional factor of e~^* (Aground -2000 g/cm )ct^n 



A 



should be included j6| in Eq. (^J. For SM interactions, this factor is very nearly unity 
and can safely be neglected but it will be important for some of the exotic models we 
will consider. 

The neutrino-nucleon cross-section in Eq. ([Q) describes both charged current (CC) 
neutrino-quark scattering and neutral current (NC) neutrino-quark scattering for which 
we adopt the cross-sections given in Ref . . The energy of the shower produced depends 
on the neutrino flavor and the type of interaction Electron neutrinos undergoing CC 
interactions produce a shower with both an electromagnetic and hadronic component: 
Esh,em = (1 — y)Eiy, -Esh,had = yE^- Muou ueutriuos undergoing CC interactions, as well 
as all neutrino flavors undergoing NC interactions, produce a hadronic shower with an 
energy, E^h^d = vEu- 

Charged current interactions of tau neutrinos are somewhat more complicated but 
more interesting. The tau lepton produced in the initial CC neutrino interaction has a 
decay length of ^ 50 m x (i?^/10^ GeV). Thus, at sufficiently high energies, the sec- 
ond hadronic shower from the tau decay will be spatially separated and be identifiable as 
a "double bang" event 0]. Lacking a full-blown simulation, we estimate that a separation 
of the two bangs by 10 km would be adequate for definitive identification by the Auger 
ground array — this requires that the primary neutrino energy exceed ~ 3 x 10^ GeV 
and that the first interaction occurs 50 km or more above the ground (which is easily 
satisfied for neutrino induced showers inclined over 80°). However, if the energy exceeds 
~ 10^° GeV, the tau lepton will hit the ground before decaying. A more careful analy- 
sis of shower profiles as seen by the Auger fluorescence detectors may allow significant 
acceptance for such events over a somewhat broader energy range. 

The scattering of electron flavor anti-neutrinos with electrons can occur efficiently via 
the resonant exchange of a W~ boson at a neutrino energy of 6.3 x 10^ GeV. Although 
we include this process in our calculations, the detector acceptance for showers at this 
energy is expected to be rather low, thus this process is of only marginal importance. 

To accurately determine the probability P[i?sh, cos6'z, X] of the Auger ground array 
observing a shower with a given energy, zenith angle and initiated at a given depth, a 
detailed detector simulation is required, which is beyond the scope of this study. To make 
a reasonable estimate we have modelled the energy dependence such that we reproduce 
the acceptances found through the simulations performed in Ref. j^. The probability 
function we arrive at is of order unity for shower energies of (9(10^^) GeV, decreases 
slowly down to energies of O(IO^) GeV, and then falls rapidly at lower energies. We treat 
hadronic and electromagnetic showers separately as in Ref. j^; in the case of a mixed 
electromagnetic-hadronic shower, we treat it as two separate showers for the purpose 
of estimating the probability of detection. In Fig. |21 we plot the spectrum of quasi- 
horizontal, deeply penetrating, neutrino induced showers as would be seen by the Auger 
ground array, for two choices of the ultra-high energy cosmic neutrino spectrum. 

C. Earth Skimming Tau Neutrinos 

A second class of neutrino events potentially observable at Auger is generated by tau 
neutrinos which interact while skimming the Earth's surface Such interactions can 
generate tau leptons which escape the Earth and produce a slightly upgoing hadronic 



shower when they decay in the atmosphere.^ This does not happen for electron neutrinos 
since the electrons produced in CC interactions are invariably absorbed in the Earth. 
For muon neutrinos, the produced muon can escape the Earth, but will not decay in the 
atmosphere since the decay length is ^ 10^ times longer than for a tau. 

When a tau lepton is generated in the Earth, it loses energy via electromagnetic 
processes at a rate per unit length of 



dEr 
dx 



-a 



(2) 



where a= 0.002 GeV cm^/g and (3 = 6x 10^^ cm^/g. At very high energies, this will often 
dramatically reduce the energy before the tau is able to decay. At moderate energies, 
this has little effect over a single decay length. In Fig. Elwe show the effect of interactions 
for tau neutrino 'Earth skimmers'. For incoming zenith angles only slightly below the 
horizon, the spectrum is not suppressed until above ~ 10^ GeV, while there is a noticeable 
pile-up near 10^ GeV. 
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FIG. 3: Effect of interactions in the Earth on cosmic tau neutrinos with a spectrum cx E, 
extending up to 10^'^ GeV. The horizontal line is the unmodified spectrum and the other lines 
are for incoming angles 0.1°, 1° and 5° degrees below the horizon. 



Tau leptons produced in CC interactions near the Earth's surface can occasionally 
escape the Earth before decaying, and thus produce a hadronic shower which is potentially 
observable by Auger. We have calculated the spectrum of tau leptons escaping the Earth's 
surface by Monte Carlo using the energy loss rate of Eq. Q . To calculate the probability 
of given tau neutrino induced shower being detected by Auger, we have used the same 
probabilities as employed in the case of quasi-horizontal showers. In addition to this 
function, however, we require that the shower be initiated at a height such that the 



^ Even when the tau lepton decays in the Earth, regeneration effects [l^ can extend the effective range, 
such that another tau lepton emerges and decay in the atmosphere. 



shower is still able to be detected (for the details of this aspect of the calculation, see 
Ref. |13|]). This is particularly important for very high energy tau leptons which can 
escape the Earth's atmosphere before decaying The Andes mountains near Auger's 
southern site are also a possible target for tau neutrinos; however, we have not included 
this effect in our calculations as the overall correction is less than 10% 3|. In Fig. E] we 
show the spectrum of Earth skimmers as would be seen by Auger. 




FIG. 4: The spectrum of Earth skimming, tau neutrino induced showers as would be seen by 
Auger. The solid and dashed lines are, respectively, for the cosmogenic neutrino flux and the 
Waxman-Bahcall flux, which would yield 1.3 and 4.8 events per year in Auger. 



III. ULTRA-HIGH ENERGY COSMIC NEUTRINO FLUXES 

Ultra-high energy neutrinos may be produced in a wide range of astrophysical sources. 
In this section, we briefly discuss some of these possibilities. 

Interactions of ultra-high energy cosmic ray protons propagating over cosmological 
distances with the cosmic microwave background generates a cosmogenic flux of neu- 
trinos nil through the decay of charged pions produced in J97 interactions jl^, which 
should also result in a suppression of the cosmic ray spectrum above the 'GZK cutoff': 
Eqzk ~ 5 X 10^° GeV. The intermediate state of the reaction P7cmb is dominated 

by the A"*" resonance, because the n decay length is smaller than the nucleon mean free 
path on the relic photons. Hence, there is roughly an equal number of 7t~^ and 7r°. Gamma 
rays, produced via 7r° decay, subsequently cascade electromagnetically on the cosmic ra- 
diation fields through e~^e~ pair production followed by inverse Compton scattering. The 
net result is a pile up of 7 rays at GeV energies, just below the threshold for further pair 
production. On the other hand, each vr"*" decays to 3 neutrinos and a positron. The e"*" 
readily loses its energy through synchrotron radiation in the cosmic magnetic fields. The 
neutrinos carry away about 3/4 of the 7r+ energy, and therefore the energy in cosmogenic 
neutrinos is about 3/4 of the one produced in 7-rays. 



The normalisation of the neutrino flux depends critically on the cosmological evolu- 
tion of the cosmic ray sources and on their proton injection spectra jl?! . 18l |. It also 
depends on the assumed spatial distribution of sources; for example, relatively local ob- 
jects, such as sources in the Virgo cluster [l^, would dominate the high energy tail of 
the neutrino spectrum. Another source of uncertainty in the cosmogenic neutrino flux 
is the energy at which there is a transition from galactic to extragalactic cosmic rays as 
inferred from a change in the spectral slope. While Fly's Eye data [2o| seem to favour a 



transition at 10^° GeV, a recent analysis of the HiRes data [21| points to a lower value of 
~ 10^ GeV. This translates into rather different proton luminosities at the sources 2^ 
and consequently different predictions for the expected flux of neutrinos jl^. A fourth 
source of uncertainty in the cosmogenic flux is the chemical composition — if ultra-high 
energy cosmic rays are heavy nuclei rather than protons the corresponding cosmogenic 



neutrino flux may be somewhat reduced 23|- Throughout this paper, we will adopt the 
cosmogenic neutrino spectrum as calculated in Ref. [l8|. 
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FIG. 5: Different possibilities for the ultra-high energy cosmic neutrino spectrum. The solid 
horizontal line (WB) corresponds to the Waxman and Bahcall bound |2a| assuming proton- 
proton interactions. The more rapidly falling solid line (WB (low crossover)) is obtained by 
the same ar gum ent but assuming a lower energy transition between galactic and extragalactic 
cosmic rays j2j| . The dotted line (Top Down) is the predicted flux in models where the highest 
energy cosmic rays arise from the decays of superheavy dark matter particles to many body 
states ji^l- Finally, the dashed line (Cosmogenic) is the spectrum of neutrinos produced in the 
intergalactic propagation of ultra-high energy protons [l^. In all cases, the curves show the 
sum of neutrinos and anti-neutrinos of all flavors. 

In addition to being produced in the propagation of ultra-high energy cosmic rays, 
neutrinos are also expected to be generated in their sources, such as gamma-ray bursts or 
active galactic nuclei . Although the details of the relationship between the cosmic ray 
spectrum and the cosmic neutrino spectrum are model dependent, some rather general 
arguments can be applied. In particular, Waxman and Bahcall j26[ have shown that for 
compact, cosmological sources, which are optically thin to proton-proton and proton- 



o 



photon interactions, an upper limit can be placed on the flux of neutrinos. We will follow 
them in adopting a neutrino spectrum arising from proton-proton collisions (with an 
inelasticity of 



El^NjdE^ < 4 X 10"^ GeVcm-2 s"^ sr" 



(3) 



summed over all flavors. After oscillations during propagation, one finds at Earth a 
nearly identical flux of the three neutrino flavors [9| with equal number of neutrinos and 
antineutrinos jl^]. If the shape of the neutrino spectrum is not an E'"^ power law, or 
if the other assumptions of the Waxman-Bahcall argument are modified, this bound can 
be exceeded [28]. For example, if their bound is evaluated under the assumption of a low 
galactic to extragalactic crossover energy (~ 4 x 10^ GeV rather than the ~ 10^'' GeV 
used by Waxman and Bahcall) a larger flux with a steeper spectrum [E'"^'^^) is obtained 
[2^ . Furthermore, sources which are optically thick such that only neutrinos can escape 
('hidden sources'), can easily exceed this bound [29l|). 

Finally, if the highest energy cosmic rays are not accelerated in distant astrophysical 
sources but are instead produced relatively locally in the galactic halo in the decays of 
supermassive dark matter particles, then significantly higher fluxes of ultra-high energy 
photons and neutrinos will also be generated S^]- This model was motivated by the 
AGASA observation that the cosmic ray spectrum continues apparently without attenu- 
ation beyond Eqik but at the same time the events are isotropically distributed on the 
sky. Such events have not yet been seen by Auger, which has moreover set a restrictive 
upper limit on the fraction of photons in the cosmic ray flux js^l- We have normalized 
the theoretical expectations for the neutrino flux from QCD and electroweak fragmenta- 
tion in heavy particle decay as calculated in Ref. [32] by matching the flux of nucleons 
observed by Auger, and checked that the photon limits are not violated. 

In Fig. El we plot the expected spectrum of ultra-high energy cosmic neutrinos in the 
models discussed above and give the corresponding event rates for Auger with standard 
QCD parton model calculations in Table H] In this calculation we have truncated the 
cosmic neutrino spectra above 10^^ GeV — this choice has only a mild effect on the 
estimated rates. 
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TABLE I: The number of neutrino induced events per year expected in Auger for various choices 
of the ultra- high energy neutrino spectrum, as shown in Figure |51 calculated using the standard 
QCD parton model cross-section. 



IV. NEUTRINO PHYSICS WITH AUGER 
A. Prospects for Cross-Section Measurements 

Deviations of the neutrino-nucleon cross-sections from the prediction of the simple 
parton model can signal new physics beyond the SM, but might alternatively be just 



n 



due to saturation effects wfiicli can substantially modify the parton density at small x (i.e. 
small energy fractions) [s^. These effects can significantly reduce the total cross-section 
at high energies, softening the power law behavior predicted by the simple 'unscreened' 
parton model toward compliance with the Froissart bound js^]- By contrast, new physics 
such as TeV-scale quantum gravity [s^ [s^ can enhance the neutrino interaction cross- 
sections.^ This has been calculated in various different frameworks, e.g., arising from 
exchange of Kaluza-Klein (KK) gravitons js^ 3^|, black hole production ji^, and TeV- 
scale string excitations 1(41,] . In this section, we will discuss the ability of Auger to measure 
deviations in the neutrino-nucleon scattering cross-section from the SM prediction, with- 
out assuming any particular interaction model. 

The event rates for quasi-horizontal and Earth-skimming neutrinos have different re- 
sponses to the inelastic cross-section j42(]. The rate of quasi-horizontal showers rises 
proportional to the cross-section (although if this exceeds ~ 10~^® cm^, attenuation of 
the neutrino flux in the upper atmosphere becomes significant By contrast, the 

rate of Earth skimming tau events is always depleted by an enhanced neutrino-nucleon 
cross-section because of absorption in the Earth. 

In order to probe deviations from the (unscreened) parton model calculation of the 
cross-section, it is necessary to note that the screening corrections affect CC and NC 
equally. To assess the experimental sensitivity to such effects, we assume a uniform 
suppression of the cross-section by a factor of 2 or 5 and show in Fig. IHlthe effect on the 
spectrum of Earth skimmers, as a function of the incoming angle. For a cross section 
reduced by a factor of 2, the total event rate of Earth skimmers is 1.4 yr~^, which is 
slightly larger than for the unscreened parton model. The reduction in cross-section due 
to screening will be energy dependent in general, but as shown in Table ITU the effect is 
mainly manifest at intermediate energies of ~ 10^ — 10^° GeV, corresponding to center-of- 
mass energies ^/s ~ 10^ — 10^ GeV; at these energies the ratio of quasi-horizontal to Earth- 
skimming events is a useful diagnostic of any suppression in the cross-section. This is in 
fact primarily because the cosmogenic neutrino flux peaks at these energies, nevertheless 
since this represents a reasonable lower limit to the expected flux, this sensitivity is 
likely to be achieved and even surpassed. A factor of 2 reduction in the cross-section 
may appear extreme, even so it is clear that Auger can probe the behavior of parton 
distribution functions (pdfs) in a kinematic region out of reach of forseeable accelerators. 
This will be particularly beneficial for callibrating different hadronic interaction models 
of air shower development, which presently differ significantly in their predictions 4^ . 

With regard to enhancements of the cross-section by new physics, in general this 
will be different for CC and for NC interactions. To assess the sensitivity of Auger, we 
consider a toy model in which only the NC cross-section is enhanced by a factor ranging 
between 3 and 100, while assuming the inelasticity to be the same as in the SM.^ In Fig. [3 
we show that this results in a suppression of Earth-skimming tau spectrum. By contrast 
the quasi-horizontal showers are enhanced, resulting in a steady increase of the ratio of 
quasi- horizontals to Earth skimmers, as the NC cross-section is increased (see Table ITTT|) . 
Clearly Auger would be sensitive to substantial increases of the NC cross-section. 

Thus both an increase and a decrease of the neutrino-nucleon cross-section from the 
naive SM value will have distinctive observational signatures. To quantitatively assess 
the sensitivity of Auger to such effects, the uncertainty in the cosmic neutrino fluxes must 



^ It is noteworthy that the neutrino-nucleon cross section can also be enhanced in some supersymmetric 

models through direct channel production of superpartner resonances js^ . 
^ This resembles the KK graviton exchange, as we discuss in the next section. 
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FIG. 6: The effect of interactions in the Earth on the tau neutrino spectrum when the (CC + 
NC) interaction cross-section is suppressed. As in Fig. |31 we adopt a spectrum oc E'"^, which 
extends to 10^^ GeV. In each frame, results are shown assuming the SM cross-section (sohd-line) 
and a cross-section smaUer by a factor of 2 (dotted-hne) and a factor of 5 (dashed-hne) . The 
three frames are for incoming angles of 1°, 2° and 5° degrees below the horizon. 



also be taken into account |4J| . Moreover, to determine the acceptances to different types 
of events, a full detector simulation is clearly required to improve over the approximate 
estimates adopted here. 
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FIG. 7: The effect of interactions in the Earth on the tau neutrino spectrum when the (NC) 
interaction cross-section is enhanced. As in Fig. 13 we adopt a spectrum cx E~'^, which extends 
to 10^2 GeV. In each frame, results are shown for the cases of the SM cross-section (sohd-Hne), a 
cross-section 10 times larger (dashed-line) and 100 times larger (dotted-line). The three frames 
are for incoming angles of 0.1°, 1° and 5° degrees below the horizon. 
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TABLE II: Variation of the rate (in yr~^) of Earth-skimming tau neutrino induced events in 
various energy intervals (in GeV), for the SM (unscreened parton) cross-section, and for a cross 
section 2 times smaller (for both CC and NC). The cosmogenic neutrino flux has been assumed. 
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TABLE III: The energy integrated rate (in yr~^) of quasi-horizontal and Earth-skimmers, as 
well as their ratio, for the SM cross-section and for different enhancements of the NC component 
alone. The cosmogenic neutrino flux has been assumed. 

B. Prospects for Flavor Ratio Measurements 

In most models of astrophysical neutrino sources, neutrinos are generated through 
the decay of charged pions: 7r+ l^^'^p. €.~^^ei^fj.^fj. or e~z/et'^t'/x, 
thus the flavor ratio at source is i/g '■ '^fj. '■ = 1/3 : 2/3 : 0. However, oscillations 
modify this ratio as neutrinos propagate to Earth. Given the observed near maximal 
mixings ji^ and the long baselines involved, the predicted flavor ratio at Earth is z/g : 
z/^ : z/t- «i 0.36 : 0.33 : 0.31 following Ref. However, cosmic (anti-)neutrinos may also 
be generated in the decay of neutrons: n — >■ p^e~i'e- In this case, the initial flavor ratio 
of z/g : z/^ : z/,- = 1 : : becomes Ue : : Ur ^ 0.56 : 0.26 : 0.18 at Earth In 
either case a measured deviation from these predictions could indicate new physics if the 
neutrino production mechanism is well understood. 

To study the sensitivity of Auger to the flavor content, we plot in Fig. |Hlthe ratio of 
quasi-horizontal showers to Earth-skimming events as the z/e flux is varied in ratio to the 
other flavors. As in Table HTH this ratio is 0.05 when the flux is equally spread among 
flavors. 



V. MODELS OF NEW PHYSICS 
A. Low Scale Quantum Gravity 

— 1/2 

Two of the most important scales in physics are the Planck scale (Mpi = ~ 
10^^ GeV) and the weak scale {Mw = Gp^^'^ ~ 300 GeV), and a long standing problem 
is explaining the hierarchy between these scales. The traditional view is to adopt Mpi as 
the fundamental scale and attempt to derive Mw through some dynamical mechanism 
(e.g. renormalization group evolution). However recently several models jl^ have 
been proposed where Mw is instead the fundamental scale of nature. In the simplest 
construction of these models, the SM fields are confined to a 3+1-dimensional 'brane- 
world' (corresponding to our observed universe), while gravity propagates in a higher 
dimensional 'bulk' space. 

If space-time is assumed to be a direct product of a 3+1-dimensional manifold and a 
flat spatial n-dimensional torus T" (of common linear size 27rrc), one obtains a definite 
representation of this picture in which the effective 4-dimensional Planck scale is related 
to the fundamental scale of gravity, Md, according to (35l | 

= 87rr;^M|^+% (4) 

where D = 4 + n. If M^ is to be not much higher than the electroweak scale, then this 
requires Vc to be large in Planck units and thus reformulates the hierarchy problem. 
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FIG. 8: The ratio of neutrino induced quasi-horizontal showers to Earth-skimming tau neutrino 
induced upgoing showers as a function of the flavor content of the cosmic neutrino flux. We 
have assumed equal numbers of muon and tau neutrinos and adopted the cosmogenic neutrino 
flux. 

For illustrative purposes in what follows we will consider only the case of flat 
extra-dimensions. The consequences of more exoteric scenarios (such as warped extra- 
dimensions [36j) have been studied in detail by various authors j47l] . 

1. Sub-Planckian Regime 

From our 4-dimensional point of view, the higher dimensional massless gravitons then 
appear as an infinite tower of KK modes, of which the lowest is the massless graviton it- 
self, while its excitations are massive. The mass-squared of each KK graviton mode reads, 
= J27=i ^11 '^1-1 where the mode numbers are integers. Note that the weakness of 
the gravitational interaction is compensated by the large number of KK modes that are 
exchanged: the coupling M^^ of the graviton vertex is cancelled exactly by the large 
multiplicity of KK excitations ~ s"/^r", so that the final product is ~ s"/^/M|,"^"' 
Here \f§ is the center-of-mass energy available for graviton-KK emission. Taking brane 
fluctuations into account, a form factor ~ e~™ is introduced at each graviton ver- 
tex ji^. This exponential suppression, which parametrizes the effects of a finite brane 
tension, provides a dynamical cutoff in the (otherwise divergent) sum over all KK contri- 
butions to a given scattering amplitude. Altogether, one may wonder whether the rapid 
growth of the cross-section with energy in neutrino-nucleon reactions mediated by spin 2 
particles carries with it observable deviations from SM predictions. 

A simple Born approximation to the elastic neutrino-parton cross section (which un- 
derlies the total neutrino-proton cross-section) leads, without modification, to ~ 

jssl . I39I I . Unmodified, this behavior by itself eventually violates unitarity. This may 
be seen either by examining the partial waves of this amplitude, or by studying the high 
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energy Regge behavior of an amplitude AR(s,t) oc s°'^^^ with spin-2 Regge pole, viz., 
intercept a{0) = 2. For the latter, the elastic cross-section is given by 

do-el \AR{s,i)\'^ £2a{0)-2 a2 

dt 

whereas the total cross-section reads 

Sm[Afi(g,0)] (o)_i 
c^tot ~ ~ s ^ ^ ~ (5j 



so that eventually, > o"tot Eikonal unitarization schemes modify this behaviour. 
Specifically, for larg e impact parameter, a single Regge pole exchange amplitude yields 
o"tot ~ ln^(s/so) |5l|, an estimate which is insensitive to the underlying theory at short dis- 
tances (UV completion). Recently, the differential cross-section for such gravity-mediated 
interaction at large distances has been calculated [s^l • Because of the large cross-sections, 
albeit with low inelasticity, there would be distinctive double and/or multiple bang 
events js^ similar to those discussed in Sec. Ill Bl [53]. For small impact parameters, 
it becomes difficult to respect partial wave unitarity as corrections to the eikonal ampli- 
tude are expected to become important. Note that graviton self interactions carry factors 
of t associated to the vertices, and thus as t increases, so does the attraction among the 
scattered particles. Eventually it is expected that gravitational collapse to a black hole 
(BH) will take place, absorbing the initial state in such a way that short distance effects 
are screened by the appearance of a horizon 4^, 5^ • 



2. Trans- Planckian Regime 



According to Thome's hoop conjecture [56[, a BH forms in a two-particle collision 
when and only when the impact parameter is smaller than the radius of a Schwarzschild 
BH of mass equal to VI = ^/xs. The total cross-section for BH production is then, 

(7) 



proportional to the area subtended by a "hoop" of radius 
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D 



M, 



D 



n+3 ' 
2 ' 



n + 2 



1 

l + n 



where F{n) is a form factor of order unity. Recent work has confirmed the validity 
of Eq. ((Tj) and evaluated the dimension- dependent constant F{n), analytically in four 
dimensions ^8] and numerically in higher dimensions js^. In the course of collapse, a 
certain amount of energy is radiated in gravitational waves by the multipole moments of 
the incoming shock waves [H^, leaving a fraction y = M^u/VS available to be emitted 
through Hawking evaporation ^]. Here, Mbh is a lower bound on the final mass of the BH 
and y/I is the center-of-mass energy of the colliding particles, taken to be partons. This 
ratio depends on the impact parameter of the collision, as well as on the dimensionality of 
space-time jH^]- Of course, this calculation is purely in the framework of classical general 
relativity, and is expected to be valid only for energies far above the fundamental Planck 
scale Md, for which curvature is small outside the horizon and strong quantum effects 
are hidden behind the horizon. Extending this formalism to center-of-mass energies close 
to Md requires a better understanding of quantum gravity. 



String theory provides the best hope for understanding the regime of strong quantum 
gravity, and in particular for computing cross-sections at energies close to the Planck 
scale j63|. In principle embedding TeV-scale gravity models in realistic string models 
might facilitate the calculation of cross-sections for BHs (and string excitations) having 
masses comparable to Mp. To be specific we will consider embedding of a 10- dimensional 
low-energy scale gravity scenario within the context of SO (32) Type I superstring theory, 
where gauge and charged SM fields can be identified with open strings localized on a 
3-brane and the gravitational sector consists of closed strings that propagate freely in 
the internal dimensions of the universe jo^]- After compactification on down to four 
dimensions, Mpi is related to the string scale, Ms, and the string coupling constant, (^s, 
by M|i = (hereafter, D = 10, i.e. n = 6). 

Subsequent to formation, the BH proceeds to decay The decay of an excited 

spinning BH state proceeds through several stages. The initial configuration looses hair 
associated with multipole moments in a balding phase by emission of classical gravita- 
tional and gauge radiation. Gauge charges inherited from the initial state partons are 
discharged by Schwinger emission. After this transient phase, the subsequent spinning 
BH evaporates by semi-classical Hawking radiation in two phases: a brief spin-down 
phase in which angular momentum is shed jH^, and a longer Schwarzschild phase. In the 
latter the emission rate per degree of particle freedom i of particles of spin s with initial 
total energy between {Q, Q + dQ) is found to be 0| 
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where Tbh = 7/(4 vrr^) is the BH temperature. 



2 7r('='-2)/2 

rp-2)/2] 



(9) 



(10) 



is the volume of a unit {d — 3)-sphere, and (Ts{Q,rs) is the absorption coefficient (a.k.a. 
the greybody factor). Recall that SM fields live on a 3-brane {d = 4), while gravitons 
inhabit the entire spacetime {d = 10). The prevalent energies of the decay quanta are 
of C(Tbh) ~ '^/^s, resulting in s-wave dominance of the final state. Indeed, as the 
total angular momentum number of the emitted field increases, as{Q,rs) rapidly gets 
suppressed |68|]. In the low energy limit, Qvg -C 1, higher-order terms are suppressed 
by a factor of 3{Qrs)~'^ for fermions and by a factor of 25{Qrs)~'^ for gauge bosons. 
For an average particle energy (Q) of C(r~^), higher partial waves also get suppressed, 
although by a smaller factor. This strongly suggests that the BH is sensitive only to 
the radial coordinate and does not make use of the extra angular modes available in the 
internal space ^]. A recent numerical study |[70i] has explicitly shown that the emission 
of scalar modes into the bulk is largely suppressed with respect to the brane emission. In 
order to contravene the argument of Emparan-Horowitz-Myers jH^ , the bulk emission of 
gravitons would need to exhibit the opposite behavior - a substantial enhancement into 
bulk modes. There is no a priori reason to suspect this qualitative difference between 
s = and s = 2, and hence no reason to support arguments [t^ favoring deviation from 
the dominance of visible decay. With this in mind, we assume the evaporation process 
to be dominated by the large number of SM brane modes. The lower bound on the mass 
radiated in the Schwarzschild phase could be somewhat reduced at large b [72] compared 
to the estimate in Ref. |62] used here. On the other hand, the effective range of b at 
which there is trapping is somewhat increase with the result that there is not any 
significant change. 



The total number of particles emitted is approximately equal to the BH entropy, 



^BH = |MBHr,. (11) 



At a given time, the rate of decrease in the BH mass is just the total power radiated 
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where Cj is the number of internal degrees of freedom of particle species i. Integration of 
Eq. dni) leads to 

Mbh = - E / ^ r(4) C(4) T^H A,, (13) 

where / = 1 (7/8) for bosons (fermions), and the greybody factor was conveniently 
written as a dimensionless constant, Fg = as{{Q),rs)/Ai, normalized to the BH surface 
area ||69| 

36 /9\2/7 2 
^4 = yvr(^-j (14) 

seen by the SM fields {Ts=i/2 ~ 0.33 and r^^i ^ 0.34 ^]). Now, since the ratio of 
degrees of freedom for gauge bosons, quarks and leptons is 29:72:18 (excluding the Higgs 
boson), from Eq. ()13p one obtains a rough estimate of the mean lifetime, 

which indicates that BHs evaporate near-instantaneously into visible quanta. 

The semi-classical description outlined above is reliable only when the energy of the 
emitted particle is small compared to the BH mass, i.e. 

Tbh <^ Mbh , or equivalently, Mbh > Mio , (16) 

because it is only under this condition that both the gravitational field of the brane and 
the back reaction of the metric during the emission process can safely be neglected (t^ . 
For BHs with initial masses well above Mio, most of the decay process can be well 
described within the semi-classical approximation. However, the condition stated in 
Eq. (jl6j) inevitably breaks down during the last stages of evaporation. At this point 
it becomes necessary to introduce quantum considerations. To this end we turn to a 
quantum statistical description of highly excited strings. 

It is well-known that the density of string states with mass between M and M + dM 
cannot increase any faster than p(M) = e^^'^^ /M, because the partition function, 

/•CO 

Zi(3) = / dMp(M) e"^^ , (17) 
Jo 

would fail to converge [t^. Indeed, the partition function converges only if the temper- 
ature is less than the Hagedorn temperature, P^^, which is expected to be ~ Mg. As P 
decreases to the transition point Pn, the heat capacity rises to infinity because the energy 
goes into the many new available modes rather than into raising the kinetic energy of 
the existing particles [t^. In the limit, the total probability diverges, indicating that the 
canonical ensemble is inadequate for the treatment of the system. However, one can still 
employ a microcanonical ensemble of a large number of similar isolated systems, each 
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with a given fixed energy E. With the center-of-mass at rest, E = M so the density of 
states is just p(M) and the entropy S = Inp(M). In this picture, equihbrium among 
systems is determined by the equahty of the temperatures, defined for each system as 



T 



dM 



M 



Equihbrium is achieved at maximum entropy when the total system heat capacity, C, is 
positive. Ordinary systems (on which our intuition is founded) have C > 0. However, 
for a gas of massive superstring excitations the heat capacity. 
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(19) 



is negative, as is the case for BHs [77]. The positivity requirement on the total specific 
heat implies that strings and BHs cannot coexist in thermal equilibrium, because any 
subsystem of this system has negative specific heat, and thus the system as a whole is 
thermodynamically unstable. This observation suggests that BHs may end their Hawk- 
ing evaporation process by making a transition to an excited string state with higher 
entropy, avoiding the singular zero- mass limit [t^. The suggestion of a string ^ BH 
transition is further strengthened by three other facts: (i) in string theory, the funda- 
mental string length should set the minimum value for the Schwarzschild radius of any 
BH jl^; (ii) Tbh ~ for ~ [so'l; (iii) there is an apparent correlation between 
the greybody factors in BH decay and the level structure of excited strings [s^l- The 
string ^ BH "correspondence principle" unifies these concepts: When the size of 
the BH horizon drops below the size of the fundamental string length 4 ^ iiQ, where 
iio is the fundamental Planck length, an adiabatic transition occurs to an excited string 
state. Subsequently, the string will slowly lose mass by radiating massless particles with 
a nearly thermal spectrum at the unchanging Hagedorn temperature j^!|.(Note that the 
probability of a BH radiating a large string, or of a large string undergoing a fluctuation 
to become a BH is negligibly small [s^].) 

The continuity of the cross-section at the correspondence point, at least parametrically 
in energy and string coupling, provides an independent supporting argument for this 
picture |63]. Specifically, in the perturbative regime, the Virasoro- Shapiro amplitude 
leads to a "string ball" (SB) production cross-section oc gls/M^. This cross-section 
saturates the unitarity bounds at g^s/M^ ~ 1 ^8^, so before matching the geometric BH 
cross-section oc r^, there is a transition region at which a ~ All in all, the rise with 

energy of the parton-parton SB/BH cross-section can be parametrized as 
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where Mio = {Sny^^Mjgl 

The inclusive production of BHs proceeds through different final states for different 
classical impact parameters b |_62]. These final states are characterized by the fraction 
y{z) of the initial parton center-of-mass energy, = ^/xs, which is trapped within the 



1 ^7 




FIG. 9: Quantitative measures of the validity of the semi-classical analysis of BH production 
for n = 6 extra dimensions, where a^min = AfBH,min/-^io- 
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FIG. 10: The cross-section for BH production in neutrino nucleon collisions, for n = 6 extra di- 
mensions, assuming Miq = 1 TeV and MBH.min = Miq. Energy losses by gravitational radiation 
have been included. The SM i/N cross-section is indicated by the dotted line. For comparison 
the typical pp cross-section is shown, as well as the cross-section required for triggering verti- 
cal and horizontal atmospheric showers. The cross-section for absorption by the Earth is also 
shown [9l| . 



horizon. Here, z = b/b^i.^, where femax = 1.3rs(v^) jH^]. With a lower cutoff MBH,mm on 
the BH mass required for the vahdity of the semi-classical description, this implies the 
joint constraint 

y{z) y/xs > MBH.min (21) 

on the parameters x and z. Because of the monotonically decreasing nature of y{z), 
Eq. (j2ip sets an upper hound z{x) on the impact parameter for fixed x. The corresponding 
parton-parton BH cross-section is cr^^^ix) = nb'^^x), where b = zb^a.x- The total BH 
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production cross-section is then j8f 

(T^^{E^, MBH,mm, Miq) = [m^ dx ^ fi{x, Q) 

) , (22) 

where i labels parton species and the fi{x,Q) are pdfs [871. '^^^ momentum scale Q 
is taken as r~^, which is a typical momentum transfer during the gravitational collapse 
process. The parameter MBH,mm plays in important role n interpreting the results derived 
below. The validity of the semi-classical calculation requires at least three criteria to 
be satisfied. First, 5*0, the initial entropy of the produced BH should be large enough 
to ensure a well-defined thermodynamic description fl^. Second, the BH lifetime tbh 
should be large compared to its inverse mass so that the black hole behaves like a well- 
defined resonance. Third, the BH mass must be large compared to the scale of the 
3-brane tension T3 so that the brane does not significantly perturb the BH metric 0]. 
Quantitative measures of these three criteria are given in Fig. |^ for n = 6, assuming 
T3 = v^87r/(27r)^ M^q for 6 toroidally-compactified dimensions ||89|]. We see that all three 
criteria are adequately satisfied for Mbh, min ~ 3Mio jl^]. The resulting uN — > BH 
production cross-section is shown in Fig. ^1 

In the perturbative string regime, i.e. MsB,min < < Mg/gs, the SB production 
cross-section is taken to be 

a^^{E^, MsB.min, Mio) = JmI^^ dx J2 fii^^ Q) o"sb(^) , (23) 

s i 

where d'^^{s) contains the Chan-Paton factors which control the projection of the initial 
state onto the string spectrum. In general, this projection is not uniquely determined 
by the low-lying particle spectrum, so there are one or more arbitrary constants. The 
analysis in the pq — >■ z/g channel illustrates this point [o^. The pg scattering, relevant 
for vN interactions at ultra-high energies, introduces additional ambiguities. In our 
calculations we adopt the estimates given in Ref. 9^1 considering the saturation limit 



and including both neutrino-quark and neutrino-gluon scattering. The resulting uN 
SB cross-section is shown in Fig. setting the Chan-Paton factors equal to 1/2. 

As can be seen in Figs. Eland although the neutrino interaction length is reduced 
below the SM value due to BH/SB production, it is still far larger than the Earth's atmo- 
spheric depth. Neutrinos therefore would produce BH/SBs with roughly equal probability 
at any point in the atmosphere. As a result, the light descendants of the BH/SB may 
initiate low-altitude, quasi-horizontal showers at rates significantly higher than SM pre- 
dictions.^ Because of this the atmosphere provides a buffer against contamination by 
hadronic showers (for which the electromagnetic component is completely attenuated at 
such large zenith angles) allowing a good characterization of BH-induced showers when 

s > 1 0,113. 



If the quasi-horizontal deep shower rate is found to be anomalously large, it can be 
ascribed either to an enhancement of the incoming neutrino flux, or to an enhancement 
in the neutrino-nucleon cross-section. However, these possibilities may be distinguished 
by focusing on events which arrive at very small angles to the horizon. An enhanced 
flux will increase both the quasi-horizontal and Earth-skimming event rates, whereas a 
large BH cross-section suppresses the latter, because the hadronic decay products of BH 



^ Additionally, neutrinos that traverse the atmosphere unscathed may produce black holes via interac- 
tions in the ice or water and be detected by neutrino telescopes |94| . 
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FIG. 11: The neutrino-nucleon cross-section including the effects of TeV scale string resonances. 
The sohd and dashed hnes correspond to models with string tension Mg = 1 and 2 TeV, 
respectively. The Standard Model cross-section is shown as a dotted line for comparison. 
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FIG. 12: The spectrum of quasi-horizontal, deeply penetrating, black hole induced showers as 
would be seen by Auger for the cosmogenic flux (left) and the Waxman-Bahcall flux (right). 
The dashed lines indicates different values of the fundamental Planck scale (from below Miq = 
10, 7, 5, 4, 3, 2, 1 TeV; in all cases MBH.min = 3Mio) while the solid line is the SM prediction. 



evaporation do not escape the Earth's crust [96j. To quantify the potential of Auger 
in discriminating BH/SB induced showers, we show separately the BH production event 
rates for quasi-horizontal and Earth skimming neutrinos in Figs. El and El The SB 
production rates are similarly given in Figs. and El and a summary of these event 
rates is provided in Tables jlVl and jVl respectively. 
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FIG. 13: The spectrum of Earth skimming, tau neutrino black hole induced showers as would 
be seen by Auger for the cosmogenic flux (left) and the Waxman-Bahcall flux (right). The 
dashed lines indicates different values of the fundamental Planck scale (from below Miq = 
1, 2, 3, 4, 5, 7, 10 TeV; in all cases MBH,min = 3Mio), while the solid line is the SM prediction. 
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FIG. 14: The spectrum of quasi-horizontal, deeply penetrating, neutrino string-ball induced 
showers as would be seen by Auger for the cosmogenic flux (left) and the Waxman-Bahcall flux 
(right). The dashed lines refer to the string scale Mg = 1 TeV (upper) and Ms = 2 TeV (lower), 
while the solid line is the SM prediction. 

B. Non-perturbative Electroweak Interactions 

The transition probability between two flat space vacua can be calculated in a 
Minkowski framework in analogy with WKB tunneling through non- vacuum fluctuations, 
or by evaluating the minimal action appropriate to a classical solution of Euclidean space 
in a given topological sector |93]. As is well known, in Yang-Mills theories the inclusion of 
massless fermions fundamentally alters the picture jo^: transitions between vacua (sep- 
arated by energy barriers whose minimum height is set by the sphaleron energy Esp 99]) 
will be totally suppressed unless accompanied by the simultaneous emission or absorption 
of all fermions coupled to the gauge field. In the Minkowski description, these fermions 
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FIG. 15: The spectrum of Earth skimming, tau neutrino string ball induced showers as would 

be seen by Auger for the cosmogenic flux (left) and the Waxman-Bahcall flux (right). The 
dashed lines refer to the string scale Mg = 1 TeV (upper) and Mg = 2 TeV (lower), while the 
solid line is the SM prediction. 
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TABLE IV: Black hole producing event rates of quasi-horizontal showers and Earth-skimming 
tau neutrino induced showers expected to be observed per year by Auger for both the cosmogenic 
neutrino flux and the Waxman-Bahcall flux. In all cases MBH,min = 3Mio. 
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TABLE V: String ball producing event rates of quasi-horizontal showers and Earth-skimming 
tau neutrino induced showers expected to be observed per year by Auger for both the cosmogenic 
neutrino flux and the Waxman-Bahcall flux. 



emerge during level-shifting in the strong 0{l/g) gauge fields interpolating between vacua 

{g = coupling constant). In the Euclidean description, the presence of a zero mode uj 
for each light fermion coupled to the gauge field will, because of the rules of Grassman 
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integration, generate a 't Hooft vertex j9 
legs, 



with all the different fermions appearing as 
uFi + h.c. , (24) 
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where Fi is a chiral fermion field. Whether these exotic processes occurwith sizeable 
rates in high energy particle collisions is a long-standing open question jlOCll | . 

At center-of-mass energies \f§ < Esp ~ irMw/cew ~ 7.5 TeV, the cross-sectio n fo r 



electroweak instanton mediated processes is known to have an exponential form [101|. 
Here, mw = 80.423 GeV is the boson mass and awimw) = 0.0338 is the SU{2) fine 
structure constant jHf. Including essential pre-exponential factors 102l |. one has, for the 



phenomenologically interesting case of fermion-fermion scattering f + f ^ all. 
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where Fhg is the "holy-grail" function [103]. By means of perturbative calculations of the 
relevant exclusive amplitudes about the instanton (/) , squaring them and summing over 
the final states, or, alternatively, by means of a perturbative calculation of the forward 
elastic scattering amplitude about the widely separated instanton anti-instanton (J/) 
pair and determining the imaginary part to get the total cross-section via the optical 
theorem, one may calculate the decisive tunneling suppressio n ex ponent Fhg, as a series 
in fractional powers of e = ^/I / {A-Kmw / aw) ^ v^/(30 TeV) [To^ . 
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Therefore, the total cross-section given in Eq. (|^^ is exponentially growing for e ^ 1. 
At e of C(l), however, the perturbative expression in Eq. (j26|l no longer apphes. In this 
energy regime, only ex trap olations of, and lower bounds on, the tunneling suppression 
exponent are available 104l |. 

Interestingly, at ~ 100 TeV, the cross-section can rise to values characteristic of 
QCD interactions. Since the electroweak instanton- induced interaction applies equally 
to all fermions, neutrinos can thus acquire hadron-like cross sections at high energies. 
Moreover, the inelasticity of the process is high. Together, these facts imply that neutrino 
interactions mediated by instantons would induce air showers in the upper atmosphere 
with characteristics similar to those of proton-induced showers |l05j . Conversely, the 
non-observation to date of deeply p enet r ating air showers constrains any sudden rise of 
the neutrino-nucleon cross-section |l06l . Il07l |. Figure ^1 shows the allowed region for 
transition from electroweak to QCD-like neutrino-nucleon cross-section, consistent with 
existing data. The dashed line indicat es th e neutrino-nucleon cross-section taken from 
Ref. [lQ§ | obtained taking o"// ^ 1 mb jl04| . As can be seen in Fig. UHl this prediction is 
marginally consistent with the region allowed by current data. For this cross- sect ion, the 
expected event rate at Auger would be 4.3 quasi-horizontal showers per year assuming 
the cosmogenic neutrino flux, and 14 quasi-horizontal showers per year assuming the 
Waxman-Bahcall neutrino flux; the rate of Earth-skimmers is 1.3 per year in both cases. 
As shown in Fig. El the suppression of Earth-skimmers due to absorption in the Earth is 
negligible. However, the rate of quasi-horizontal showers is increased by about 2 orders of 
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FIG. 16: The allowed 90%, 95%, and 99% CL regions for interpolation between the electroweak 
and QCD-like neutrino-nucleon cross-section consistent with existing data. Also shown with a 
dashed line is t he p r edic ted enhancement of the cross-section by electroweak sphalerons. For 
details see Ref. [lo3,Eoi- 



magnitude, and such events would be concentrated in a small energy range, as indicated in 
Fig. ^1 This would provide a clean signal for electroweak instanton-induced interactions. 
Thus, if no deeply developing showers are observed, tighter constraints can be placed on 
this model, and more generally on any sudden rise in the neutrino-nucleon cross-section. 
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FIG. 17: Suppression of Earth skimming events due to electroweak sphalerons as would be seen 
by Auger for the cosmogenic neutrino flux (left), and for the Waxman-Bahcall flux (right). The 
solid line is the SM prediction. 



C. Neutrino Decay 

Neutrinos are known to be sufficiently light that they are stable against tree-level 
electroweak decays. Moreover, decays of the form z/j uj'-f ot u ^ vvV are severely con- 
strained by experiment However, some models of lepton number violation postulate 
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FIG. 18: The spectrum of quasi-horizontal showers mediated by electroweak sphalerons as 
would be seen by Auger for the cosmogenic neutrino flux (left), and for the Waxman-Bahcall 
flux (right). The solid line is the SM prediction. 



the existence of a massless Goldstone boson, the Majoron, X. Consequently, decays such 
as Ui — >■ VjX or Vi — >■ VjX, are then possible, where denote mass eigenstates 109]. 
Presently such possibilities are only weakly constrained by Solar neutrino data, which set 
the bound r/m ^ 10~^ s/eV [llO|. However, by studying cosmic neutrinos which have 
travelled over far longer baselines Auger can be more sensitive to their instability by a 
factor of ~ 10^ — 10^, if an effective flavor ratio measurement can be made. Because of 
the extremely large energies pro bed by A uger, it will be complementary in this regard to 
the IceCube neutrino telescope 111, lid . 



The ratio of flavors observed in the cosmic neutrino spectrum depends on whether any 
species of neutrinos have decayed and on the decay channel. In the simple situation where 
all heavy neutrino species decay into the lightest mass eigenstate (or into non-interacting 
states, such as a sterile neutrino), we would expect to observe at Earth the flavor ratio 

0.. : (t>u, ■■ (pu. = cos^ 9q : ^ sin^ 9^ : ^ sin^ ^0^6:1:1, (27) 

where 6q is the solar neutrino mixing angle and we have assumed the normal hierarchy 
as well as Ues = 0. This result is independent of the flavor ratio at source. However, for 
the case of an inverted hierarchy, the predicted flavor ratio at Earth is 

0.. : K ■■ = : U% : [7,^3 ^ : 1 : 1, (28) 

where [/q- is the neutrino mixing matrix and we have taken the atmospheric mixing 
angle to be maximal. These results are in striking contrast to the expectation for stable 
neutrinos discussed earlier in Sec. IIVBI 

These two cases represent the most extreme deviations from the usual phenomenology 
and are robust in that they do not depend on the flavor composition at source. A 
variety of other (more baroque) possibilities have been considered, e.g. only the heaviest 
neutrino eigenstate decays, but the p redic ted flavor ratios after propagation then depend 



on the assumed flavor ratio at source |llll |. In Table IVIl we list some of these possibilities 
assuming the usual mass hierarchy and source flavor ratios as for pion decay. 

We cannot measure the flavor ratios directly at Auger. However, as discussed earlier. 
Earth-skimming events are generated uniquely by tau neutrinos, while quasi-horizontal 



showers can be generated by all neutrino flavors. Furthermore, because of maximal mixing 
of z/^ and Vr we expect their fluxes to be always comparable. Therefore, by combining 
these two measurements Auger can potentially determine the flavor ratios of ultra-high 
energy neutrinos. 



Decaying Mass Eigenstates 


Decay Products 




1^3, ^3 


Irrelevant 


6:1:1 


1^3 


Invisible 


2:1:1 


1^3 




1.4-1.6:1:1 


1^3 




2.4-2.8:1:1 


^3 


50%iyi, 50% 1/2 


2:1:1 



TABLE VI: The neutrin o fla.v or ratios predicted for a variety of neutrino decay models with 
decay mode as indicated |lll | . 



Of course this requires a substantial event rate. As shown in Table H] the standard 
cosmogenic neutrino flux is expected to generate only about 0.7 quasi- horizontal shower 
events over 10 years. This is certainly insufficient for making the precision measure- 
ments needed to identify the effects of neutrino decay. For the nominal Waxman-Bahcall 
flux. Auger is expected to detect about 2.2 quasi-horizontal events and about 48 Earth- 
skimming events in 10 yr. How ever, if the cosmic ray galactic-extragalactic transition 
happens at around 10^ GeV |ll3j , then the required proton luminosity in the extragalactic 
sources increases signiflcantly. Then Auger would detect as many as 21 quasi-horizontal 
and 350 Earth-skimming events over 10 years. This corresponds to a 2a measurement of 
their ratio of 0.06 ± 0.026, which would exclude anomalous flavor composition with a z/g 
content greater than : (j)^^ : (j)^^ ~ 2.5 : 1 : 1. 

Other possibilities for altering neutrino flavor ratios have been explored 114l . Uli . If 
Lorentz invariance is violated through modiflcation of the usual dispersion relation for 
neutrinos by non-renormalizable operators induced by quantum gravity effects, then the 
fraction of tau neutrinos may be suppressed 115l |. Auger would then observed the ratio 
of Earth-skimmers to quasi-horizontal events to decrease from about 20 to close to zero. 



VI. CONCLUSIONS 

Our knowledge of fundamental interactions has largely been limited to the energies 
up to which collider experiments have been able to probe. The Tevatron, currently 
operating at Fermilab, produces collisions with a center-of-mass energy slightly below 
2 TeV, while the Large Hadron Collider, under construction at CERN, will reach 14 
TeV. By contrast, a typical neutrino observed at the Pierre Auger Observatory will have 
an energy of O{10^) GeV, corresponding to a neutrino-nucleon center-of-mass energy 
exceeding 40 TeV. Although the number of collisions which will be observed (i.e. the 
beam luminosity) is far below that of collider experiments. Auger and other experiments 
sensitive to ultra-high energy cosmic neutrinos have in principle the ability to provide 
unique information on new physics beyond the reach of any planned accelerator. 

In addition to this advantage, cosmic neutrinos have traveled over very great distances 
before reaching Earth, thus their detection also constitutes an extremely long-baseline 
oscillation experiment. Instead of being limited to phenomena which occur over minuscule 
fractions of a second, cosmic neutrinos provide an exceptional window into phenomena 
only evident over cosmological scales of length or time. 



The Pierre Auger Observatory is capable of detecting two primary classes of neutrino 
induced events — quasi-horizontal, deeply penetrating showers and (slightly) upgoing 
showers induced by Earth-skimming tau neutrinos. Used separately, the rates of such 
events are of limited use in probing new physics; since the spectrum of cosmic neutrinos 
is currently unknown, an event rate cannot by itself be used to determine the neutrino- 
nucleon interaction cross-section. However by combining these two classes of neutrino- 
induced events, it becomes possible to make a crude cross-section measurement. As 
this cross-section is increased (decreased) , the rate of quasi-horizontal showers increases 
(decreases) accordingly, while by contrast, the rate of shghtly upgoing showers is reduced 
(enhanced) since Earth- skimming tau neutrinos become absorbed more (less). Thus the 
ratio of quasi-horizontal, deeply penetrating showers to slightly upgoing showers provides 
a check of the behaviour of the neutrino-nucleon cross-section at ultra-high energies. The 
details of such a measurement, of course, depend on the energy dependence of such 
interactions, as well as their inelasticity and other characteristics. 

These two types of neutrino-induced events also provide the opportunity to constrain 
the ratios of flavors present in the ultra-high energy cosmic neutrino spectrum. If these 
neutrinos are generated through the decay of charged pions (as they are in most models), 
they will reach Earth in nearly equal quantities of each flavor after oscillations are taken 
into account. A larger than expected rate of quasi-horizontal, deeply penetrating showers, 
in comparison to the slightly upgoing shower rate, would thus indicate a suppression of 
the tau neutrino component in the ultra-high energy cosmic neutrino spectrum, due, for 
example, to neutrino decay. 

In this study, we have considered several specific models in which either the neutrino- 
nucleon cross-section, or the ratio of cosmic neutrino flavors, deviates substantially from 
the expectation of the perturbative Standard Model. We have studied enhancements in 
the neutrino-nucleon cross-section in models with low-scale gravity, variously described 
as due to the exchange of Kaluza-Klein gravitons, the production of microscopic black 
holes, and/or string resonances. We have also considered increases in the cross-section due 
to non-perturbative Standard Model electroweak instanton induced processes, which in 
contrast do not lead to a decrease in the inelasticity. Regarding flavor ratio measurements, 
we have discussed several models of decaying neutrinos. 

It is difficult to precisely delineate the reach of these techniques as this depends on 
the unknown flux of cosmic neutrinos at the energies to which Auger is sensitive. We 
have considered both the "guaranteed" cosmogenic flux which sets a lower bound and 
the Waxman-Bahcall flux which sets an upper bound. Further observations of ultra-high 
energy cosmic rays by Auger itself will help to pin down the expected neutrino flux. 

Over the next few years, the Pierre Auger Observatory may well identify the world's 
flrst ultra- high energy neutrino event. We have attempted to illustrate the exciting new 
possibilities for probing new physics that will be opened up by such a detection. 
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